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摘 E 本 文 首先 从 内 积 的 角度 给 出 了 移动 最 小 二 乘 逼 近 法 的 新 的 推导 方法 ， 然 后 对 Lancaster 等 提出 的 
移动 最 小 二 乘 播 值 法 进行 了 重新 推导 ， 取 在 插值 节点 奇异 的 权 函 数 ， 并 对 基 函 数 进行 部 分 正 交 
化 ， 建 立 了 改进 的 移动 最 小 二 乘 插值 法 ， 并 证 明了 其 形 函 数 的 插值 性 质 。 本 文 提出 的 改进 移动 最 
小 二 乘 插值 法 的 公式 比 Lancaster 的 公式 更 为 简单 ， 并 可 提高 形 函 数 的 计算 效率 。 本 文 为 工程 问 
题 的 插值 型 无 网 格 方法 提供 了 建立 形 函 数 的 基本 方法 。 
关键 词 : 移动 最 小 二 乘法 ， 移 动 最 小 二 乘 插值 法 ， 形 函数 ， 权 函数 ， 紧 支 域 
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1 引言 


无 网 格 方法 是 继 有 限 元 法 之 后 的 又 一 类 重要 的 科学 和 工程 分 析 的 数值 方法 。 移 动 最 小 二 
FREIA (Moving least-squares approximation) 是 无 网 格 方法 构造 形 函 数 的 最 为 重要 的 方法 之 
一 所。 由 于 移动 最 小 二 乘 训 近 法 的 形 函 数 不 具 有 Kronecker Delta 函数 的 性 质 ， 因 而 基于 该 方 
法 建立 的 无 网 格 方法 需要 借助 罚 函 数 法 或 Lagrange HF iA MMA IRA, SBE NSS 
形式 方程 较为 复杂 ， 降 低 了 无 网 格 方法 的 计算 效率 。 所 以 ， 研 究 具 有 播 值 性 质 的 移动 最 小 二 乘 
插值 法 (Interpolating moving least-squares method) 是 非常 重要 的 。1981 F, Lancaster 等 在 进 
行 曲 面 拟 含 时 推广 了 标准 的 最 小 二 乘 形式 ， 提 出 了 移动 最 小 二 乘 逼近 法 和 移动 最 小 二 乘 插值 
法 阅 。 移 动 最 小 二 乘 通 近 法 的 形 函 数 比 移动 最 小 二 乘 插值 法 的 形 函 数 简单 得 多 ， 以 至 于 后 来 发 
展 的 基于 移动 最 小 二 乘 概念 的 无 网 格 方法 基本 上 都 是 基于 移动 最 小 二 乘 逼近 法 。 

Mukherjee 等 对 移动 最 小 二 乘法 进行 了 改进 ， 以 方便 无 单元 Galerkin 方法 处 理 边界 条 件 里 ， 
陈 美 娟 和 程 玉民 等 提出 了 改进 的 移动 最 小 二 乘法 外 ， 程 玉民 等 提出 了 复 变量 移动 最 小 二 乘 
法 加 及 其 相应 的 复 变 量 无 网 格 方法 6- 引 。 本 文 首先 从 内 积 的 角度 给 出 了 移动 最 小 二 乘 逼 近 法 的 
新 的 推导 方法 ， 然 后 对 Lancaster 等 提出 的 移动 最 小 二 乘 插值 法 进行 了 重新 推导 ， 取 在 插值 节 
点 奇异 的 权 函 数 ， 并 对 基 函 数 进行 部 分 正 交 化 ， 建 立 了 改进 的 移动 最 小 二 乘 播 值 法 ， 并 证 明了 
其 形 函 数 的 插值 性 质 。 本 文 提 出 的 改进 的 移动 最 小 二 乘 插值 法 的 公式 比 Lancaster 的 公式 更 为 
简单 ， 可 提高 形 函数 的 计算 效率 。 
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2 ”移动 最 小 二 乘法 
EER RRO 上 的 函数 w(x)， 已 知 其 在 域内 N 个 节点 xl, x2,… ,xm 的 函数 值 ， 取 函数 9 
uh(x) = > pi(x) -a;(x) = p™(x)-a(x), x e Q, (1) 
i=l 
为 函数 u(x) 的 逼近 函数 。 这 里 Pi(x) 是 基 函 数 ， 痉 是 基 函 数 的 个 数 ，ui(x) 是 相应 的 系数 。 基 
函数 的 通常 形式 为 


线性 基 : 
p” =(1,21), 对 于 一 维 区 域 . (2) 
p7 = (1,z1,z2)， 对 于 二 维 区 域 . (3) 

二 次 基 
pz =(1,21,22), 对 于 一 维 区 域 . (4) 
p7 = (1, 21,22,22, 2122,22), 对 于 二 维 区 域 . (5) 


对 应 于 (1) 式 的 整体 通 近 ， 在 点 x 的 邻 域内 的 局 部 逼近 定义 为 
u(x, x) = 》 pi(X) .ai(x) = p” (x) - a(x). (6) 
i=l 


系数 ai(x) (i = 1,2,… ,m) 可 由 局 部 近似 的 加 权 最 小 二 乘 拟 合 得 到 ， 即 对 每 一 点 x, a(x) 的 选 
择 总 是 使 下 列 定义 的 泛 函 取 极 小 值 。 定 义 


J = X w(x- x) [u (x, xr) — u(x1)}? 


= Yuk- xr) | narn s(x) uber] (7) 


I=1 i=1 


其 中 w(x — xr) 是 具有 紧 支 集 特 性 的 权 函 数 ，xr (7 = 1,2,---,n) 为 点 x 的 紧 支 域内 的 节点 。 
(7) 式 可 用 矩阵 形式 表示 为 


J = (pa — u)T W(x)(pa — u), (8) 
其 中 
u= (uug sUn)? = (ua) ula) +++ u(x)) (9) 
pi(x1) p2(X1) + pm(X1) 
= ae een _ ii = (P1; P2,°** Pm); (10) 


pi(xn) pe(Xn) …: Pm(Xn) 





第 6 期 任 红 萍 等 : 改进 的 移动 最 小 二 乘 插值 法 研究 1023 
w(x — x1) 0 toe 0 
Wee we Mee s (11) 
0 0 “++ w(x — Xn) 
a = a(x) = (a:(x), as (x), ,an(x)) (12) 
对 于 函数 f(x)， 在 点 x 我 们 记 
E = (f(x1), fo),  f(%n))”, (13) 
定义 函数 f(x), g(x) 在 x 点 的 内 积 为 
(F, 9)x = (f, 8)x = £7 W(x)g = > w(x — xr)f(xr)g(x7), (14) 
I=1 
相应 地 ， 定 义 在 x 点 的 范 数 为 ， 
ll = [(£, f)<] °: (15) 
H (8) 式 得 
J = |ipa — ull? = |laipi + a2p2 + +--+ @mPm 一 ullx”, (16) 
QW u 4E28 lH] span(pi, p2,… , Pm) 上 的 投影 向 量 u(1) 对 应 的 
J = jw -ul (17) 
BIAS). BCR u IEA OMEN 
E me} (18) 
其 中 
ul) € span(pi,D2,:'' Pm), (19) 
u?) £ span(pı, P2;:: 5 Pm); (20) 
所 以 
(Pi, u)x = (Pi， a), + (Pis u) = (Pi; u)),., (21) 
即 
(Pi, u)x = (Pp; u),. = (P, Zup) (22) 
j=1 
写成 矩阵 形式 ， 即 为 
(Pl,pl)x (P1,P2)x … (P1,Pm)x ay(x) (pi,u)x 
(p2, pi)x (pe, Po) x Sees (Pe, Pm)x ao(x) (p2, u)x 
. : . . i = í (23) 


(Pm; pi )x (Pm; P2)x WY (Pm, Pm)x am (x) (Pm, u)x 
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BJ) 
A(x)a(x) = B(x)u, (24) 
这 里 
A(x) = PT W(x)P, (25) 
B(x) = PT W(x). (26) 
解 方程 组 (24) 可 得 
a(x) = A !(x)B(x)u, (27) 


把 式 (27) RAR (1) 即 得 到 逼近 函数 的 表达 式 


u" (x) = B(x)u = > *;(x)ur, (28) 
I=1 
ZE P(x) 称 为 形 函 数 ， 且 为 
B(x) = (8, (x), Bo(x), +++ , Bn(x)) = p” (x)A `! (x)B(x). (29) 


当 m = 工时 ， 唯 一 的 基 函 数 pi;(x) = 1, HEIN 
w(x — xr)ur 





u(x) = E. (30) 
> w(x — xr) 
YA HAR BR 
v(x 一 xi) = ale Se) i=1,2,---,n, (31) 
> w(x — xy) 
id 
v(x) = (u(x — x1), v(x 一 x2) : ° ,0(K—Xn))”, (32) 
us(x) = Sox —x,)u(xy) = v? (x)u. (33) 


3 改进 的 移动 最 小 二 乘 插值 法 


移动 最 小 二 乘法 得 到 的 和 逼近 函数 是 不 经 过 插值 节点 的 ， 若 要 逼近 函数 经 过 所 有 插值 
节点 ， 就 必须 使 权 函 数 w(x -xi) 在 xi 点 处 是 奇异 的 ， 但 这 就 导致 方程 组 (24) 的 系数 矩 
阵 A(x) = P7W(ojJP 在 x 一 xi 时 趋 于 co， 为 了 解决 这 个 问题 ， 我 们 对 基 函 数 进行 如 下 正 
在 空间 span(pi, p2,… Pm) Po ARM pi (x) 三 1 在 x 点 单位 化 为 


p= Pre orni (34) 
| Pa P [ 37 w(x- xr)] ° 


I=1 
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再 将 P2,P3,°°° ,Pm 与 bw 在 x 点 正 交 化 为 


2 w(x — xr)pi (xr) 





BL) (x) = pi(x) — (pi, bE?) bY? = pi(x) 


5 w(x ~ x7) 
Tt 
一 pi (x) iay Y v(x — XI)pi(X1), t= 2,3, , Tn, (35) ` 
I=1 
Bp 
b® (x) = pi(x) — p? (x), (36) 
i 
WP (x1) Poa) 1 0 Ox) 
(2) (3) ee (m) 
Gua p E a (37) 
BP (xn) WS (xn) Q OL” (xn) 


现在 我 们 对 新 的 基 BD (x), BP (x), --- DL? (x) 使 用 移动 最 小 二 乘法 ， 注 意 到 基 中 后 面 m — 
1 个 与 第 一 个 之 间 是 正 交 的 ， 以 及 第 一 个 的 范 数 是 1， 则 逼近 函数 为 


uh (x) = (u, BY?) <b (x) + Da iC 0), (38) 
= 
即 
u(x) = v (x)u + bT (x)a(x), (39) 
其 中 
B(x) = (P (x), bE (x), + BL 0), (40) 
a(x) = (a1 (x), a2(x); -+ , am-1(x))”. (41) 
由 移动 最 小 二 乘法 可 得 
a(x) = AZI(x)Bx(x)u， (42) 
这 里 
Ax(x) = CTW (x)Cx, (43) 
Bx(x) = CT W(x). (44) 
将 式 (42) 代入 式 (39) 可 得 
u®(x) = vT(x)u + bT(x)A-1(x)B, (x)u. (45) 
sÑ (45) 可 以 写成 
ub (x) = B(x)u = > B17(x)ur, (46) 


I=1 
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更 (x) = v” (x) + b’ (x)Ax!(x)Bx(x). (47) 

将 式 (42) 与 Lancaster 得 到 的 公式 
a(x) = AX! (x)Bx(x)(u — (v’ (x)u)p,), (48) 


相 比 ， 式 (42) 更 简洁 。 


4 ”改进 的 移动 最 小 二 乘 插值 法 的 插值 特性 


引 理 1 车 取 权 函数 为 i 





w(x — Xi) 


~ |x — x, |’ 


其 中 a ATER, ¿= 1,2,---,n, v(x- xi) 由 式 (31) 给 出 ， 则 有 


n 
1) 0<v(x 一 xi)< 1; 2) SS u(x —xr) = 1; 3) lim v(x — x;) = 6; 
T=1 XX; 


4) lim w(x — xk)u(x — xj) = w(xj — xx), k £ j; 
X— Xk 








5) lim f%(x)=f(x:); 6) lim bẸ®(x) = 0; lim b(x) = 0. 
区 一 Xi X—x; x—x; 
证 明 1), 2) BR. PUES), 47 =i 
Ü | 
lim v(x — x;) = lim I ee = — = 1, 
x—x; eC Fe 
当 了 天 1 时 
1 
lim %(x 一 xi) = lim poe 0. 
XX; sey ER s ee ee oe l: 
Ix—xil® Ilx—x;|® |x—xnl® 
对 于 4)， 当 k 关 j 时 
lim w(x — xk)u(x — x;) 
X— Xk 
| > e] | 1 
= lim xk T <s dle T = —ə T = w(Xk = xj). 
AI eget aaa aka eee |%;k 一 x; |° 
对 于 5)， 因 为 


n 


f°(x) = >》 u(x ~ xr)f (xi), 


I=1 


故 有 


要 证 6)， 由 式 (36) 并 注意 到 (5) 可 得 


lim f*(x) = lim $ v(x- x1) f(r) = lim v(x- x) fa) = f(x). 
,全 : 


lim 02) (x) = py(xi) — lim pš(x) = p; (xi) — p; (x) = 0. 
x—x; XX 
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7) 由 式 (40) 及 (6) 可 得 。 证 毕 

因为 权 函 数 在 插值 点 的 奇异 性 ， 我 们 先 计算 经 过 正 交 处 理 后 的 基 函 数 之 间 的 内 积 以 及 它们 
和 任意 函数 的 内 积 当 自 变量 趋 于 插值 点 时 的 极限 情况 。 

引 理 2 `x — xk hf, (b) u). M (bY, bi) G = 2,3,--- ,m) 的 极限 存在 ， 且 分 别 为 


1) im (bŠ, u)x = lm) — pi(xx)]u (x; — xx) [u(x7) — u(x,)]; 


2) im (be, bx))x = jag POD — pi(xk)]ju (xr — xx)|[pi(xr) — pi(xk)]- 


故 当 x 一 xk 时 ，Ax(x), B,(x) 的 极限 存在 ， 从 而 a(x) 的 极限 存在 。 
证 明 1) 首先 由 内 积 定义 知 


(bQ, u)x = > bD (xr) (x g; xr)u(xr), 





I=1 
由 引 理工 的 结论 可 知 
| 5 w(x — x;)pi(;) 
bY (xer)w(x — xul) = (pier) =a w(x — xz)u(xz) 
> w(x — x;) 


II 


u(xr)o(x =x) S> o(x—x;)(pi(xr) — pilx;)), 


j=l GAL 


所 以 


m n 
im (OK) Wx = dim XY u(xro(x~xr) SS w(x —x;)(pi(xr) — pi(x;)) 
T=1 j=1,j4I 


|. 


Jim ys u(xr)u(x — xy)w(x — xx) (pi(xr) — pi(xk)) 
I=1,IZ#k 
+u(xk) >` w(x; — xk)(pi (xk) — pi (x;)) 
j=1j#4l 


n 


> (pi(xr) 一 Pi(xh))w(xr — xr) (u(x) — u(xg)). 


I=1,I¢k 


lI 


2) 由 内 积 定义 知 


n 
PP, BO) = J P Cen) w(x — x7)0P (xr), 
I=1 
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由 (1) 及 引 理 1 的 结论 有 


lim (b® pb), = lim Dee- xr) ii w(x — xj) (pi(xr) — pi(x;)) 


XX, X—xk 
j=l gAl 
n 
= (p(x/)-pi(xe))u(x;— xe) (pilxr) — p)). 
I=1,IZk 


定理 Jim uh (x) = u(x;) = ui. 
证 明 HK (39) 以 及 引 理 2 得 


lim u"(x) = Jim u" (x) + lim b7(x)- lim a(x) = u(x;). 
一 Xi x—x; XX; 


其 -一 Xi 


5 ”结论 


本 文 首先 从 内 积 的 概念 给 出 了 移动 最 小 二 乘 逼 近 法 的 新 的 推导 方法 ， 然 后 对 Lancaster 等 提 
出 的 移动 最 小 二 乘 插值 法 进行 了 改进 ， 取 在 插值 节点 奇异 的 权 函数 ， 并 对 基 函 数 进行 部 分 正 交 
化 ， 建 立 了 改进 的 移动 最 小 二 乘 插值 法 ， 并 证 明了 其 形 函 数 的 插值 性 质 。 本 文 提出 的 改进 的 
移动 最 小 二 乘 插值 法 的 公式 化 Lancaster 的 公式 更 为 简单 ， 可 提高 形 函 数 的 计算 效率 。 基 于 本 
文 方法 建立 的 无 网 格 方法 可 以 方便 地 直接 施加 本 质 边界 条 件 ， 减 少 以 往 大 多 数 借助 罚 函 数 法 
或 Lagrange 乘 子 法 施加 本 质 边 界 条 件 时 求 弱 形 式 较为 复杂 的 问题 ， 从 而 提高 无 网 格 方法 方程 
的 计算 效率 。 
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Researches on the Improved Interpolating Moving 
Least-squares Method 


REN Hong-ping!, CHENG Yu-min?, ZHANG Wu? 


(1- School of Applied Science, Taiyuan University of Science and Technology, Taiyuan 030024; 
2- Shanghai Institute of Applied Mathematics and Mechanics, Shanghai University, Shanghai 200072; 
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Abstract: In this paper, based on the concept of an inner product, we propose a new method to obtain 
the moving least-squares approximation. The interpolating moving least-squares method proposed by 
Lancaster is thus improved. By taking a singular weight function in the interpolating points and 
orthogonalizing some of basis functions, an improved interpolating moving least-squares method is 
formulated and the property of the corresponding shape function is proved. The new method has a 
simpler formula and higher computing efficiency compared with the Lancaster’s method. The proposed 
method can be used to form the shape function of the meshless method for engineering problems. 

Keywords: moving least-squares approximation; interpolating moving least-squares method; shape 


function; weight function; compact support domain 
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